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ABSTRACT 

We present deep imaging at 6.7/^ni and 15/im from the CAM instrument 
on the Infrared Space Observatory {ISO), centred on the Hubble Deep 
Field {HDF). These are the deepest integrations pubhshed to date at these 
wavelengths in any region of sky. We discuss the observation strategy 
and the data reduction. The observed source density appears to approach 
the CAM confusion hmit at 15/zm, and fluctuations in the 6.7/zm sky 
background may be identifiable with similar spatial fluctuations in the 
HDF galaxy counts. ISO appears to be detecting comparable field galaxy 
populations to the HDF, and our data yields strong evidence that future 
IR missions (such as SIRTF, FIRST and WIRE) as well as SCUBA and 
millimetre arrays will easily detect field galaxies out to comparably high 
redshifts. 

Key words: galaxies: formation ~ infrared: galaxies - surveys 



1 INTRODUCTION 

The Infrared Space Observatory {ISO, Kessler et al. 
1996) offers an improvement of orders of magnitude in 
sensitivity over the IRAS satellite, at least at shorter 
wavelengths (< 20/im). ISO is expected both to sam- 
ple the intermediate redshift {z ~ 0.5 — 1) star forming 
galaxy population at lower, less extreme luminosities, 
and to detect more strongly star forming galaxies to 
far higher redshifts. We used the CAM instrument 



on ISO (Cesarsky et al. 1996) in Director's discre- 
tionary time to observe the Hubble Deep Field (HDF, 
Williams et al. 1996), resulting in the deepest surveys 
to date at 6.7/xm and 15/im. The obvious advantage 
of this field is the extensive multi-wavelength follow- 
ups either published or underway {e.g. Dickinson et 
al. 1997 in preparation, Fomalont et al. 1996, Cowie 
et al. 1996) as well as unparallelled deep multicolour 
optical morphologies from the Hubble Space Telescope 
{HST). 
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Parameter 


LW-2 6.^^^m 


LW-3 IS^im 


Pixel field of view 


3" 


6" 


M, N steps 


8, 8 


8, 8 


M, N step size 


5", 5" 


9", 9" 


Tint 


10 sec 


5 sec 


■'Vstab 


80 


100 




10 


20 



Table 1. Summary of our CAMOl Astronomical Observa- 
tion Templates (parameters defining our M X N rasters). 
An exposure of Tint is made A^obs times at each of the 
M X N raster positions. Each raster was proceeded by 
Astab readouts at the first raster position to stabilise the 
CAM detector. 



The HDF galaxy population is strikingly domi- 
nated by blue objects which may be comparable to 
local giant Hll regions, and which are consistent with 
significant star formation. This is re-enforced by their 
largely disturbed structures. In such an interpretation 
a substantial fraction of the luminosity from young 
massive stars is absorbed by dust and re-radiated in 
the mid- and far-infrared. Clearly, imaging in the mid- 
infrared samples the spectral energy distributions of 
such galaxies much closer to this significant, and per- 
haps dominant, contribution to the bolometric power 
output. 

In this paper we present our CAM images of 
the HDF, discuss the data reduction steps and 
make crude comparisons of the sky fluctuations with 
smoothed HDF images. Subsequent papers will ad- 
dress the source identification algorithm (paper II, 
Goldschmidt et al. 1997); the source counts and com- 
parison with models (paper HI, Oliver et al. 1997), 
including a more sophisticated treatment of the con- 
fusion noise; the associations with HDF galaxies (pa- 
per IV, Mann et al. 1997); and the spectral energy 
distributions of our sources and implications for star 
formation history (paper V, Rowan-Robinson et al. 
1997). 

The maps presented here are available at 



http : //artemis .ph. ic. ac .uk/hdf or from the au- 



thors. 

2 DATA ACQUISITION AND ANALYSIS 
2.1 Observation Strategy 

The observing strategy was designed to reach the pre- 
dicted confusion limit of CAM, using the source count 
models of Pearson & Rowan- Robinson (1996). The 
observations were performed in two bands, 6.7/um and 
15^m (the LW-2 and LW-3 filters respectively), to ob- 
tain (albeit limited) colour information. These bands 
are the widest available in CAM so yield the deepest 
possible integrations. Typically more than one HDF 
galaxy falls within the Airy disk, even at 6.7/im, so 



sub-pixel offsets were needed to maximise the avail- 
able spatial information. 

We made three 8x8 rasters in microscanning 
mode (CAMOl) analagous to the "dithering" in the 
HDF, each centred on one of the HDF Wide Field 
(WF) frames, at both 6.7/im and 15/im, using in total 
~ 44.9ks of time including overheads. Our choice of 3" 
(6") pixel sizes at 6.7/im (15/im) and raster step sizes 
of 5" (9") yields optimal fiat fielding accuracy and 
spatial resolution. At 6.7/.im this gives a 96" field of 
view, well-matched the HDF WF frames. The space- 
craft jitter observed in-flight is ±0.5" {2a limits, half- 
cone), much smaller than our choices of either pixel 
size or step size. Cosmic ray transients (discussed be- 
low) ruled out readout integration times longer than 
10 seconds, and the signal-to-noise vs. time predic- 
tions from the CAM simulator implied diminishing 
returns for more than 10 (20) readouts per raster po- 
sition at 6.7/im (15^m). The number of stabilisation 
readouts prior to the rasters at each wavelength is 
appropriate for 5cr sources. The Astronomical Obser- 
vation Template (AOT) parameters are summarised 
in Table 1; for a more detailed discussion of these pa- 
rameters see the CAM Observers Manual^ Resulting 
noise levels in each of the fields are listed in table 2. 
The pixel scale of ISOPHOT made longer wavelength 
observations impracticable. 

The edited raw data FITS files (supplied by ESA) 
were processed using the CIA (CAM Interactive Anal- 
ysis April 1996 version) IDL package with the excep- 
tion of the deglitching, the construction of the fiat 
field and the mosaicing of the rasters, as discussed 
below. 



2.2 Dark subtraction, deglitching and flat 
fielding 

The default dark frame in the April 1996 CIA version 
was subtracted from the data. 

Cosmic ray events were easily identified in the 
readout histories of each pixel as > 4(t rises followed 
(one or two readouts later) by > 4cr falls. A similar 
algorithm was used to find readout troughs. These 
events were masked out in the mosaicing discussed 
below. 

However, a minority of cosmic rays appear to 
cause transients in subsequent readouts with roughly 
exponential decays (see figure 1) persisting over a few 
readouts. These glitch transients are in general diffi- 
cult to model. No attempt was made to identify and 
remove them; instead, they are effectively removed by 
median filtering in the mosaicing below. 

Use of the ESA-supplied flat field gave very un- 
satisfactory results. Instead, we created our own sky 



The 



CAM 



manual 



available 



from 



tittp: / /isowww. estec . esa.nl/inanuals/iso_cajn/ 
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Band 


A 


Width 


Area 


Date 


Jmodel 


-fiso 


la/ 


Normalisation 




(^tm) 


(Aim) 




Observed 






beam 


Shift/ Add 


Drizzle 


LW-2 


6.75 


3.375 


WF-2 


28/6/1996 


73.4 


77.3 


7.78 


1.03431 


1.03390 








WF-3 


28/6/1996 




79.9 


8.59 


1.01630 


1.01492 








WF-4 


28/6/1996 




80.0 


8.43 


1.00000 


1.00000 


LW-3 


14.5 


4.833 


WF-2 


1/7/1996 


508 


401 


36.8 


1.00000 


1.00308 








WF-3 


1/7/1996 




402 


27.4 


1.00255 


1.00000 








WF-4 


26/6/1996 




396 


46.1 


1.01306 


1.01454 



Table 2. Summary of Calibration and Sky Backgrounds. Intensities quoted are in /ijy per arcsec^. The sky backgrounds 
in each of the fields (/iso) S'l'e compared against the model (/model) discussed in the text. Noise estimates assume a 
circular beam size of 6" for LW-2 and 12" for LW-3. 



flat by noting that each detector pixel samples 64 dif- 
ferent sky positions during the raster. For each detec- 
tor pixel, wo examined the histogram of (unmasked) 
readouts and fitted Gaussians to find a mean value. 
To eliminate both sources and glitch transients, > 5a 
outliers from the mean were eliminated and the fit 
was iterated. 



2.3 Shift-and-add maps 

The rasters wore mosaicked together using two com- 
peting algorithms: drizzling (described below), and 
variants on shift-and-add. The latter are expected to 
have higher signal-to-noise but at the expense of spar 
tial resolution. 

In order to eliminate possible long timescale de- 
tector sensitivity drifts, the rasters were renormalised 
by the following method. First, we calculated the me- 
dian readout in each of the 32 x 32 CAM pixels for 
each raster at each wavelength. Second, we found the 
moan in the central 11 x 11 of these images. This 
determined the relative renormalisations, which are 
listed in Table 2. Comparison of the source positions 
in each individual raster revealed a systematic offset 
in the WF-4 Ib/j-m frame, probably due to a random 
offset in the lens positioning. For consistency with the 
drizzled mosaics we applied the same offset to this 
raster as discussed below for the drizzling. 

In the shift-and-added frames, we began by defin- 
ing an image with pixel size one-sixth that of the 
CAM detector pixels (ie, one-sixth of 3" at 6.7/im, 
and of 6" at 15/^m), encompassing the area surveyed 
by all three rasters. Each position in this fine-gridded 
image may have been observed several times by the 
CAM detector array, so we compiled a list of such 
CAM pixel readouts for each fine-gridded image po- 
sition. The data at a given sky position could then 
(for example) be median filtered or mean averaged to 
produce a final image. 

The glitch transients discussed above make a 
large contribution to the noise in final mosaics made 
by simple mean averages of readouts. An obvious al- 
ternative is median filtering, though this has a signal- 
to- noise penalty (about \/2). The glitch transients 



have timescales of order or less than the duration of 
a pointing, so a possible compromise is to mean aver- 
age over pointings, and median filter the means. This 
may also have the advantage of identifying affected 
readouts more efficiently. 

Several mosaics were therefore created from the 
readout arrays: (a) using the median readout at each 
position, to eliminate glitches; (b) using the mean 
readout at each position, with the ±5cr outliers elimi- 
nated and iterated to eliminate glitches; (c) using the 
mean readout within each raster position, followed 
by the median of these means; (d) using a clipped, 
iterated mean of pointing means. We found the best 
6.7/ttm map to be the simple median filtered image, 
but at 15/ttm the optimal map is the median-of-means, 
presumably reflecting a greater sensitivity to glitch 
transients in the final mosaic. 



2.4 Drizzled maps 

The drizzled images were produced with the same 
code (Pruchter & Hook, 1996) used to produce the 
optical HST images. Briefly, instead of superimpos- 
ing overlapping pixels, the drizzle algorithm allows 
the user to shrink the input pixel sizes (the "foot- 
print" ) before superimposing. At one limit the drizzle 
algorithm is equivalent to interlacing; at the other 
it is similar to shift-and-add above. Since the ISO 
HDF images were taken with fractional pixel spac- 
ings between them maps with increased resolution can 
be constructed. The footprint was chosen to produce 
roughly the highest resolution image without leaving 
gaps in the output image. 

Glitch transients were first removed from the 
input images using the following median filtering 
scheme, noting that at 6.7/^m the raster step offsets 
correspond to l| pixels and at 15/xm the offsets are 
I5 pixels. The 64 individual pointings were grouped 
into sets whose members have integer pixel offsets (9 
groups for LW-2 and 4 for LW-3), and for each of these 
groups the medians of the data at each pixel position 
(allowing for the shifts) was put into an image and 
the variance was computed. These 9 images at 6.7/tim 
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and 4 images at served as input to the drizzling 

algorithm. 

The satellite astrometry (RA, Dec and spacecraft 
roll angle) from the individual pointings was com- 
bined into these (9+4) data sets by taking the average 
of those pointings making up each group, after shift- 
ing them by the nominal raster offsets. This astrom- 
etry was used to define the tangent-plane projection 
used in the final maps. 

The ISO images were drizzled assuming there 
were no geometrical distortions in the images so that 
rows and columns were always assumed to be paral- 
lel and equally spaced. The 6.7/im images were sub- 
sampled to 1/3 of the input pixel size, resulting in 
output pixels of 1". The input pixel footprint was set 
to 0.65 and the pixel weights were taken from the 
reciprocals of the variances. The 15/im images were 
sub-sampled to 1/2 of the input pixel size, resulting 
in output pixels of 3". The footprint was again set to 
0.65 and the same weighting scheme was used. 

The drizzle code preserves flux by sharing the in- 
put levels amongst the output pixels. To return the 
output images to intensity units the 15/jm images 
were multiplied by 4 and the 6.7/xm images by 9. 
The input weights were the reciprocals of the vari- 
ances, and to convert the output weights from the 
drizzle code back into a variance a multiplying factor 
of 1/footprint^ (2.367) had to be applied. These out- 
put variances are not totally independent as each has 
contributions from neighbouring pixels. 

The 15Atm (LW-3) image for field WF-4 ap- 
pears to have incorrect astrometry in that when the 
appropriate offset was applied the sources did not 
line up exactly with the sources in field WF-2 or 
WF-3. An empirical offset was calculated using the 
correl_images routine from the IDL astrolib librar y 
( ittp : // idlastro .gsf c .nasa.gov/homepage . html ) 
to cross- correlate the images and from visual check- 
ing. The resultant shift was 6.4" in the detector x- 
direction and —0.7" in the y-direction. 

To ensure the final mosaics were not affected by 
global variations in the background the medians from 
the drizzled images of the individual pointings were 
compared and a multiplying factor (listed in Table 2) 
was applied to all the frames making up a pointing 
to bring them to the same level. The final drizzled 
images were then recomputed. 

2.5 Flux calibration 

To estimate the sky background, we constructed his- 
tograms of pixel counts for the central third of the un- 
normalised drizzled maps of each individual field. We 
fitted Gaussians to each of these histograms, and the 
resulting sky levels are listed in Table 2. The relative 
sky levels are in excellent agreement with the rela- 
tive normalisations calculated above. A slightly more 
sophisticated procedure was adopted for the noise es- 
timates. We selected a grid of positions in the central 



~ 1/3 of each drizzled mosaic, and placed a circular 
aperture of the Airy disk size on each position. We 
could then estimate the noise level on the scale of the 
Airy disk from the histogram of counts enclosed by 
these apertures. 

Flux calibration assumes the standard conver- 
sion in the CAM handbook. In Table 2 we compare 
our background measurements with a zodiacal back- 
ground model. This model linearly interpolates be- 
tween entries in Table 8 of the CAM Observer's Man- 
ual, but assumes a 275K blackbody zodiacal spectrum 
(Hauser et al. 1984) rather than the alternative in the 
manual. The solar elongation angle at the time of our 
observations was approximately 71.3°. The cirrus con- 
tribution at both wavelengths is expected to be < 2 
per cent of the zodiacal background. 

The sky calibration appears on the whole to be 
quite satisfactory. Checks on the calibrations using 
the sky, such as these, do not of course check the lin- 
earity of the response to sources, which may be very 
difficult to address with the ISO HDF data. It may 
nevertheless be possible to estimate this by examin- 
ing the average temporal profiles of detected sources. 
Further discussion on the sensitivity to point sources 
is contained in paper II (Goldschmidt et al. 1997). 



3 RESULTS 

3.1 Comparison of maps 

Figures 2 and 3 show the results of the drizzling al- 
gorithm. The correspondence between the drizzling 
and the shift-and-add mosaics (not shown) is excel- 
lent suggesting that neither algorithm is introducing 
artifacts. Figures 4 and 5 show the coverage maps at 
the two wavelengths. Note the poorer coverage at the 
edges, which causes the poorer signal-to-noise at the 
corresponding edges of the mosaics. 

3.2 Source confusion 

Several sources are clearly seen at 15/^m, and the 
source identification will be discussed in a later pa- 
per. However, we note here that a visual inspection 
clearly shows that we are approaching the confusion 
limit of one source per 40 beam sizes. 

One surprising feature in the 6.7/im maps is the 
apparent sky fluctuations at around 3% of the back- 
ground, which may be caused by source confusion, or 
by glitch transients which have been smeared out by 
the mosaicing process. We tested this by taking ra- 
tios of each of the 6.7/im shift-and-add mosaics. One 
apparent (albeit dubious) source was less prominent 
with increasing median filtering, so is a good candi- 
date for such a smeared-out transient. Apart from this 
the ratios of the 6.7/im shift-and-add mosaics are flat 
to ~ 1% accuracy, and do not have structures corre- 
sponding even roughly to the observed 6.7/im maps. 
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This suggests that the fluctuations are not an arte- 
fact of transients or flat-field errors, but are due to 
marginally detected sources. 

A more convincing demonstration was made by 
shuffling the CAM pixels at random and making mo- 
saics. The resulting frames should be equally suscep- 
tible to smeared-out transients, but genuine sources 
and background structure should be dispersed evenly 
over the mosaic. These randomised frames show much 
less structure than the real frames {e.g., figure 6), 
again suggesting that this structure is not an instru- 
mental artefact. 

Note however that this does not test for the pres- 
ence of correlated pixel-to-pixel fiat field variations. 
If such transient effects are present, they would be 
indistinguishable from genuine sky structure in the 
mosaicked images, and would be extremely difficult 
to remove in the data reduction process without also 
removing genuine structure. (It is nevertheless not 
clear how such correlated fiuctuations would arise.) 
One approach to testing this is to rotate or refiect the 
arrangement of the CAM detector pixels, but repeat 
the same mosaicing as before. Any sky structure will 
be dispersed over the image (though not necessarily 
evenly), but the resulting images would be equally 
sensitive to structure from correlated pixel-to-pixel 
fluctuations. These images should therefore show less 
structure than the correctly mosaicked images. Fig- 
ure 6 shows an example of this test, and the frame 
does indeed show less structure than the correctly mo- 
saicked frame. However, it is not clear how much of 
the remaining structure is due to unevenly dispersed, 
genuine sky structure, and how much is due to cor- 
related pixel noise. These potential artefacts are still 
under investigation, but in the meantime it should be 
noted that some of the structures may potentially be 
due to instrumental artefacts. 

Finally, a comparison of suitably smoothed driz- 
zled and HST F814W images, discussed in the next 
section, shows remarkable correspondence. 



3.3 Field distortion 

To compare the apparent low level sky fiuctuations 
with the underlying HDF galaxy distribution, we con- 
volved the drizzled ISO mosaics with an empirical 
point spread function (Oliver et al. 1997, Goldschmidt 
et al. 1997). In figures 7 and 8 we show a greyscale 
reproduction of the F814W HDF and flanking fleldsjj]. 
Contours of our smoothed drizzled ISO mosaics are 
overlaid. Two important points should be noted: flrst, 
the apparent low level structures do indeed match the 
inhomogenous galaxy distribution in the HST F814W 



t The flanking field image (unfortunately without astrom- 
etry) is obtainable from 

http: / / www.stsci.edu/ftp / science /hdf/project /flanking. html 



frame, at least at Ibfim; second, this match is not per- 
fect, suggesting an unknown astrometric error. 

The astrometric errors do not appear to be ran- 
dom, but rather systematic and varying continuously 
over the field. This is unlikely to be due to telescope 
pointing and lens positioning errors, since these would 
produce a non-varying systematic offset (recall that 
the individual rasters were registered with respect to 
each other). A likely explanation for the slight astro- 
metric errors is therefore field distortion in CAM, i.e. 
the CAM detector array is not quite square. In-fiight 
calibration of this distortion by the CAM consortium 
is currently underway; preliminary results indeed sug- 
gest that the total field distortion is about one pixel 
over the whole array, and that the distortion is con- 
tinuous. Once this distortion is quantified we expect 
to release appropriately revised mosaics. In the mean- 
time, further papers in this series will conservatively 
treat the astrometry as subject to random errors of 
order the Airy disk size. 

At 6.7/im the stucture, if present, is close to 
the noise; nevertheless, there are clearly several 
marginally detected sources at the positions of bright 
galaxies in the HDF F814W image: for example, sev- 
eral bright galaxies in field 2 (top left HDF frame) lie 
on or close to contour peaks. Note that the signal-to- 
noise decreases sharply towards the edges, resulting in 
many features without apparent I-band counterparts. 
Finally, recall that correlated pixel-to-pixel fiat field 
variations (if present) would mimic real sky structure. 
Nevertheless, at least some of the stucture appears to 
be marginally detected or unresolved sources. 



4 DISCUSSION 

Visual inspection of the maps at 6.7^m and 15/^m 
clearly shows many sources and these will be dis- 
cussed further in the following papers. At these wave- 
lengths the radiation is increasingly dominated by the 
re-processing of starlight by dust and it is to be ex- 
pected that the detected galaxies have above average 
star-formation rates. Our data demonstrates that ISO 
can detect galaxy populations comparable to those in 
faint optical surveys, and may also have implications 
for proposed confusion limited surveys at around our 
wavelengths, such as that from the WIRE satellite. 
However, we defer a more extensive discussion of the 
confusion limit to paper HI (Oliver et al. 1997). 

These results are also promising for surveys at 
other wavelengths. Mobasher et al. (1996) combined 
their HDF photometric redshift database with a star- 
burst spectral energy distribution model and the HDF 
galaxy counts, to obtain predictions for cumulative 
number counts of HDF galaxies at wavelengths long- 
ward of 6.7/im. Their expectation of ~ tens of sources 
at ISO wavelengths at our approximate limiting flux 
densities per beam is clearly seen to be broadly cor- 
rect. A detailed comparison with source count models 



© 0000 RAS, MNRAS 000, 000-000 



6 Serjeant, Eaton, Oliver, et al. 



is included in paper III (Oliver et al. 1997); however, 
we note here that Mobasher et al. also used the same 
models to predict ~ a hundred sources in the HDF 
at both 60/im and 0.8mm to ~ 10 — lOO/xJy levels. 
Our results are thus clearly encouraging for SCUBA 
and the proposed large milimetre arrays, as well as for 
planned future infrared space missions such as FIRST 
and SIRTF. 

Further information on the ISO HDF project 
can be found on the ISO HDF WWW pages: 
tittp : //artemis .ph. ic. ac.uk/hdf/ 
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Figure 1: Example glitches. The data shown is a dark subtracted pixel history with the (logarithmic) y-axis in the 
instrumental aiialoguc-to-digital units per second (ADU/s) and x-axis in readouts. Note the clear glitches as well as 
the bright glitch at readout ~ 580 with a transient. 

Figure 2: Drizzled mosaic at 6.7/im. North is approximately 10 degrees right of vertical and 
east is to the left. 

Figure 3: Drizzled mosaic at 15/im. 

Figure 4: Coverage map at 6.7/im. Note the poor coverage at the edges. The scale converts the grey levels to the 
relative coverage. 

Figure 5: Coverage map at 15/im. As with the 6.7/im counterpart, the coverage is poor at the edges. 

Figure 6: Comparison of the 6.7/im drizzled image (left) with a similar (shift-and-add) mosaic made after randomising 
the CAM detector pixel positions (right), and a drizzled mosaic obtained after rotating the CAM detector array 

through 90° (centre). As argued in the text, the drizzling and shift-and-add methods (not compared here) arc found 
to be in excellent agreement. Note the lack of structure in the randomised frame, and the less prominent structure 
in the rotatcd-detector frame. As discussed in the text, the latter attempts to disperse any true structure on the 
sky over the image, while being equally sensitive to correlated pixel noise. 

Figure 7a; The H3T F814W image overlayed with contours of the smoothed drizzled 6.7/(m image discussed in the 
text. 

Figure 7b: The smoothed drizzled 6.7/tm image discussed in the text, with contours overlaid. 

Figure 8: The HST F814W image overlayed with contours of the smoothed drizzled 15/im image discussed in the 
text. 
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